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Recent experiments indicate that the spin-density waves (SDWs) in (TMTTF) 2 Br, (TMTSF) 2 PF 6 
and a-(BEDT-TTF)2MHg(SCN)4 are highly unconventional and coexist with charge-density waves 
(CDWs). We present a microscopic theory of this unusual CDW-SDW coexistence. A complete 
understanding requires the explicit inclusion of strong Coulomb interactions, lattice discreteness, 
the anisotropic two-dimensional nature of the lattice, and the correct bandfilling within the starting 
Hamiltonian. 



PACS indices: 71.30.+h, 71.45.Lr, 75.30.Fv, 74.70.Kn 

Organic charge transfer solids (CTS) exhibit a bewil- 
dering array of exotic phases at low temperature, includ- 
ing superconductivity (SC), SDW, CDW, and bond-order 
wave (BOW) and (related) spin-Peierls (SP) phases. 
Among the intriguing parallels between organic and high 
temperature superconductors is the proximity of the SC 
phase to the SDW [0J|]. Given the enormous focus on 
understanding the "normal" state of high temperature 
superconductors, it seems clear that an analogous effort 
should be made to understand the "normal" state of or- 
ganics that exhibit SC. The SC-SDW proximity then sug- 
gests that it is crucial to understand the precise nature 
of the SDW and the mechanism of its formation, as well 
as the roles played by dimensionality effects and strong 
Coulomb interactions |l|,^|. Very recent experiments on 
several materials indicate that the SDWs in the organics 
are highly unconventional in nature , im- 

posing strict requirements of any theory of the SDWs. 

In this Letter, we present a unified theory of SDW 
formation in a large class of 2:1 cationic CTS. While we 
believe that our theory is general, we focus on the specific 
materials (TMTSF) 2 PF 6 , (TMTTF) 2 Br and (BEDT- 
TTF) 2 MHg(SCN) 4 (M = K, Rb, Tl), in which the un- 
conventional nature of the SDW has been demonstrated 
recently. In the past, experiments by several groups have 
established that the low temperature insulating phases 
in (TMTSF) 2 PF 6 fo| and (TMTTF) 2 Br 0] are SDWs. 
Surprisingly, recent X-ray scattering experiments ||^|] 
have revealed features associated with CDW in both 
materials even for T < T SDW . In (TMTTF) 2 Br, the 
experiment finds signatures of even a 4kj? lattice dis- 
placive instability (kp = Fermi wavevector), along with 
the more usual 2kp charge instability within the SDW 
phase The status of the experiments in (BEDT- 

TTF) 2 MHg(SCN) 4 is very similar: while early mag- 
netic susceptibility fll2[ and /isr [ [T3| measurements pro- 
vided clear evidence for a SDW, measurement of angle- 
dependent magnetoresistance oscillations fH|,@,|| have 
engendered the view that the insulating phase here is 
a "mysterious" state that is a "SDW accompanied by 
CDW" or a "CDW accompanied by SDW" @ . The 



authors of a recent 13 C-NMR study have concluded that 
the insulating state here is not a SDW at all, but is a 
CDW ||. The authors, however, also state that they 
have "no idea which kind of CDW reconciles the suscep- 
tibility anisotropy ... and other magnetic properties." 

Coexistence of CDW-SDW of the same periodicity is 
outside the scope of standard theories of density waves 
in the organic CTS |lj,|l^], which emphasize the nest- 
ing associated with the quasi-one dimensional (quasi-ld) 
Fermi surface of these materials. Although CDW, SDW 
and SP instabilities are possible within these theories, 
these instablities occur in nonoverlapping parameter re- 
gions, making coexistence impossible. The microscopic 
theory we present here explains the puzzling CDW-SDW 
coexistence in a natural fashion. The attractive features 
of our theory include: (i) our Hamiltonian is the standard 
strongly correlated model for quasi-ld organic conductors 
p6[ , and although the manifestations are novel, no ex- 
otic interactions are necessary to generate them; (ii) the 
theory can explain the differences between structurally 
similar quasi-ld systems that exhibit the SP phase and 
quasi-2d systems that exhibit the SDW; and finally, (iii) 
the theory clarifies the limitations of single-particle nest- 
ing concepts Jl^] and of the theoretical bias that SDW 
exists only at or near the 1/2-filled band. 

We posit that the Hamiltonian appropriate for the ma- 
terials considered here is the quasi-2d extended Hubbard 
model, 
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In the above, j is a site index while M is a chain in- 
dex, Bj, ktLt M,a = [Cj,L,cr C k,M,cr + fl.C.], Aj <M = {Uj,M ~ 

u j+i,Ai), where Uj t M is the displacement of the molec- 
ular site from its equlibrium position, and v^m repre- 
sents an intramolecular vibration. The total Hamiltonian 
describes coupled chains, with on-site Coulomb interac- 
tion U, intrachain nearest-neighbor Coulomb interaction 
V, and intra- and interchain nearest-neighbor hoppings 
tj = t — a(uj t M ~ Wj+i,a/) and t±. For simplicity, we 
assume a rectangular lattice. We are interested in the 
realistic parameter regime t± ~ 0.1, V ~ 2\t\, U > 4\t\. 

A critical implicit parameter in the Hamiltonian of 
Eq. (1) is the bandfilling. Charge-transfer from the 
cations to the inorganic anions in the 2:1 CTS of interest 
leads to one hole per two organic molecules, i.e., a 1/4- 
filled band of holes. This is the bandfilling we consider. 
Based on the weak dimerization along the stack axis (ob- 
served even above the metal-insulator transition temper- 
ature Tmi P]||)) it is commonly argued that 2:1 cationic 
CTS can be modeled as effective quasi-ld 1/2-filled band 
systems jl4|,[l5| with Fermi surface nesting. We show that 
although some aspects of the physics of the strongly cor- 
related 1 / 4- filled band can be understood within the weak 
coupling effective 1/2-filled theory, others simply cannot. 
In particular, since the coexistence of CDW-SDW with 
the same periodicity is impossible within the 1/2-filled 
band 0, the recent observations ||,|,||J|0,|,§ clearly 
preclude this scenario as a consistent description of the 
normal state of these CTS. We postpone further discus- 
sion of this issue until later, when we show that both from 
theoretical and experimental perspectives, the weak high 
temperature dimerization cannot lead to a modulation of 

An intuitive understanding of the broken symmetry 
coexistence within Eq. (1) can be obtained in the Id 
limit, t± = 0. Since long-range SDW does not oc- 
cur here, the relevant order parameters are the site 
charge density and the bond-order IJjj+i,M,M,<7-' Pe- 
riodic modulations of the charge density lead to the 
CDW, and of the bond-order to the BOW. Previ- 
ous work has established |lq] that the periodic lat- 
tice distortion arising from the BOW has the form 
uj = uo[r2cos(2kpj — 82) + r^cos{Akpj — 64)], where r-i 
and r4 are the relative weights of the 2k F and Ak F 
components, and 82 and 84 are the corresponding 
phase angles [|l8|. In contrast, the CDW can have 
either the 2k F or the Ak F modulation but not both, 
so rij = 1/2 + no(cosQj — (j)),, where Q — 2kp or Ak F 
jl^fni. For comparison with what follows, we sketch 
in Figs. 1(a) and (b) the familiar BOW and SDW con- 
figurations for the 1/2-filled band. A crucial feature of 
non-l/2-filled commensurate bands is the symbiotic co- 
existence between the BOW and CDW jS[L9|; as shown 
in these references, this coexistence also enables one to 
understand the general case from the f3 = limit, which 
we henceforth adopt. Importantly, the 1/4-filled band is 



unique among the non-l/2-filled bands in that Coulomb 
interactions can drive BOW-CDWs that are different 
from that driven by electron-phonon interaction. For in- 
stance, as the Hubbard U is increased from zero, the 
phase angles (82, 4>) of the 2kp BOW-CDW switch from 
(0,7r/4) to (7r/4,0) Q. Wc show the uncorrected and 
correlated 2k F BOW-CDWs in Figs 1(c) and (d), re- 
spectively. For nonzero V < V c (where V c — 2\t\ for 
U — > 00, and is larger for finite U) the absolute ground 
state acquires a Akp BOW character [r^ 7^ 0, 84 = 0, see 
Fig. 1(e)], but the CDW continues to have periodicity 
2k p |T^,|2C]]. Our numerical results will establish that the 
BOW-CDWs of Figs. 1(d) and (e) can coexist with the 
particular SDW shown in Fig. 1(f), consistent with the 
experimental data. 

We begin our numerical analysis in the Id limit. Since 
a true long-range SDW cannot occur here, we incorpo- 
rate an additional (external field- like) term Hsdw — 
— J2j ^[nj^cos(2k F j) + n.j^cos{2kpj + ip)] and consider 
the ground state of H + Hsdw- Hsdw imposes a SDW 
in the 1/2-filled band for ip = tt and the SDW of Fig. 1(f) 
in the 1/4-filled band for -0 = ir/2, with the amplitude of 
the SDW increasing with e. We examine the BOW-SDW 
coexistence using exact finite size calculations. Specifi- 
cally, we calculate the exact electronic ground state en- 
ergies E(auo, e) of finite periodic rings as functions of e, 
where auo is a rigid bond modulation parameter. The 
quantity AE(auo,e) = E(omq = 0,e) — E(auQ 7^ 0,e) 
is a direct measure of the energy gained on bond dis- 
tortion. For the 1/4-filled band, it is also necessary to 
specify ^/r^. While we have confirmed that our results 
are valid for many different r2 and n, we show in the fol- 
lowing the specific cases of r2 = (2kp BOW-CDW only) 
and r2 = r4 (both 2kp and 4k p, contributing equally). 

In Fig. 2 we show the behavior of AE(auo, e) for a 1/2- 
filled band of 10 sites and a 1/4-filled band of 16 sites, for 
U = 6, V = 1. Because of the larger energy gain the 1/2- 
filled band, the nonzero values of auo are 0.05 (in units 
of \t\) in the 1/2-filled case and 0.1 in the 1/4-filled case, 
respectively. AE(ccuq, e) decreases rapidly with e in the 
1/2-filled band, in agreement with the known result that 
the SDW and the BOW are mutually exclusive here [O . 
In contrast, we find that AE(auo 7 e) increases with e in 
the 1/4-filled band for both the cases studied, indicat- 
ing a cooperative interaction between the BOW and the 
SDW. Since the BOWs we have studied coexist with the 
correlated 2kp CDW, the ground state of H + Hsdw is 
an admixture of (2k F + 4A: F )-BOW, 2k F CDW and 2k F 
SDW for the 1/4-filled band. We have further determined 
that the bond orders for nonzero e show a modulation of 
the type shown in Fig. 1(e) for the 1/4-filled band even 
for zero lattice distortion (uq = 0), indicating a tendency 
for spontaneous BOW distortion in the presence of the 
SDW. 

For Hsdw = a true SDW can occur only for t± ^ 0. 
We have therefore performed calculations of spin-spin 
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correlations, site charge densities and bond-orders in the 
ground state of H alone for coupled chains. These cal- 
culations were done using the constrained path quantum 
Monte Carlo (CPMC) approach @. In CPMC a con- 
straining wavefunction (in this case the free-electron solu- 
tion) is used to approximate nodal boundaries and avoid 
the sign problem present in 2D calculations. We expect 
that sign problems are less severe for the 1/4 filled band 
and for the small t± we have used. All calculations were 
checked against exact results for a 4 x 4 lattice. The 
CPMC calculations are for four coupled chains of length 
twelve sites each, with the same values of U, V, auo as in 
Id and t± =0.1. The 12x4 lattice is taken to be periodic 
along both directions. We incorporate a phase difference 
of 7r between the BOWs on neighboring chains, based on 
calculations (a) in the noninteracting limit, and (b) for 
the 4x4 lattice in the interacting cases, that indicate 
that this particular phase difference gives the lowest to- 
tal energy ||. A SDW coexisting with the BOW-CDW 
requires now antifcrromagnetic interchain spin-spin cor- 
relations. In Fig. 3 we show the spin-spin correlations be- 
tween consecutive sites 2,3,4,5 on the first chain and sites 
1-12 on the second chain for the case r 2 = only. As seen 
from the figure, (a) The 2Uf bond distortion leads to anti- 
ferromagnetic interchain spin-spin correlations, (b) there 
is a simultaneous intrachain antiferromagnetic spin-spin 
correlation, - the spin densities on sites 3 and 4 are op- 
posite to those on sites 2 and 5, and (c) the magnitude 
of the interchain spin-spin correlation for a given site on 
the second chain does not simply decrease with the sep- 
aration from the site on the first chain, but is also de- 
termined by the charge density on the particular site on 
the second chain. This is a signature of long range SDW 
within the distorted lattice. With our choice of t±, the 
CPMC technique does not allow us to obtain sufficiently 
accurate spin-spin correlations between sites two chains 
apart. However, with a slightly larger t± (0.2), these 
more distant interchain spin-spin correlations are also in 
agreement with antiferromagnetic interchain correlations 
. Finally, as in Id, we have also performed the calcula- 
tions for r 2 = and find that the interaction between the 
BOW-CDW and the SDW remains cooperative, decreas- 
ing (increasing) auo decreases (increases) the strength of 
the interchain antiferromagnetic correlations p2| . 

In Fig. 4 we sketch the ground state broken symmetry 
that emerges from the CPMC results: two adjacent sites 
with unequal charge but parallel spins are surrounded by 
other such pairs with opposite spins. Viewing the pairs 
of sites as a single effective sites, this appears similar to 
the SDW of the effective 1/2-filled band scenario 
but there is the critical distinction that there are differ- 
ent charge and spin densities on the individual molecules 
within the pairs, and this internal structure will show up 
in experiments. In their X-ray diffraction experiment, 
Pouget and Ravy find strong evidence only for the 2k p 
CDW in the so-called SDW phase of (TMTSF) 2 PF 6 , and 



perhaps both the CDW as well as a 2k p or Akp BOW in 
the SDW phase of (TMTTF) 2 Br §[§. Our results indi- 
cate that the CDW-SDW necessarily drives a modulation 
of the tj. Whether or not an observable lattice distortion 
accompanies this merely depends on the strength of a. 

We now revisit the issue of the bandfilling. The param- 
eters in Eq.(l) should be derived from the overall Hamil- 
tonian describing both the organic cations and the inor- 
ganic anions. The crystal structures of (TMTTF) 2 X and 
(TMTSF) 2 X indicate that the anions face the "stronger 
bonds" between consecutive organic molecules ||. The 
anions therefore introduce a nonnegligible "anionic po- 
tential," the leading term in which takes the form 
u [ c j m <j c j+t,M,<j + h.c], where v (=1) is the number of 
(extra) electrons on the anion. This interaction modifies 
the intrachain tj , decreasing the effective hopping integral 
between the particular pair of organic molecules, equal- 
izing consecutive tj, and stabilizing a metallic phase. 
Very similar conclusions have been reached by Fritsch 
and Ducasse ^jj, who actually demonstrated the de- 
crease in the hopping integral in question by taking the 
anion potential into account. Note that equal effective 
hopping integrals at T > Tmi is essential for explaining 
the high temperature metallic conductivity if Coulomb 
correlations are strong: any dimerization of the hopping 
integrals leading to an effective 1/2 filled upper subband 
should have, in the presence of a nonzero U, opened a 
Mott-Hubbard gap, leading to an insulator. 

In conclusion, we have shown that the "normal" 
ground state of the 2:1 cationic CTS is a complex mix- 
ture of density waves. In particular, we have discov- 
ered a cooperative interaction among the CDW-BOW 
and the SDW, which emerges naturally when lattice dis- 
creteness, Coulomb interactions and actual bandfilling 
are taken properly into account. Three final comments 
are in order. First, in the appropriate small t± regime, 
our theory correctly describes the existence of highly Id 
systems (such as (TMTTF) 2 X, where X^Br) which ex- 
hibit only SP/BOW and CDW coexistence (i.e., no long- 
range SDW) (22). Second, although lack of space pre- 
cludes detailed consideration of the important magnetic 
field-induced SDW phenomena ^4|, the discussion sur- 
rounding Fig. 4 suggests that key features of the pre- 
vious approaches fl5|| remain true within the 1/4-filled, 
strongly correlated framework. In addition, however, in- 
teresting effects due to unequal charge and spin densities 
on the paired sites may emerge. Third, since SC ap- 
pears in these CTS only upon the melting of the SDW, it 
seems that any theory of organic SC should take into con- 
sideration the important roles of bandfilling and strong 
Coulomb interactions that are established by the present 
work. 
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FIG. 1. Schematics of the Id (a) 1/2-filled BOW, (b) 
1/2-filled SDW, (c) 1/4-filled uncorrelated 2k F BOW-CDW, 
(d) 1/4-filled correlated 2k F BOW-CDW (the 4k F BOW is 
the same as (a)), (e) 1/4-filled BOW-CDW that is a super- 
position of 4k F (Fig. 1(a)) and 2k F (Fig. 1(d)) periodicities, 
and that occurs for nonzero V < V c , (f) 1/4-filled 2k F SDW 
investigated here, and (g) the 4k F site-diagonal CDW that oc- 
curs only for V > V c . The double (dotted) bonds are strong 
(weak), the single bond is of intermediate strength, and the 
double dotted bond in (e) is a weak bond that is stronger 
than the single dotted bond. The heights of the vertical bars 
(arrows) on sites give their relative charge (spin) densities. 
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FIG. 2. AE(au ,e) versus e for (a) the Id 1/2-filled band 
(solid line), (b) the Id 1/4-filled band with the bond distortion 
of Fig. 1(d) (long dashed line), and (c) the Id 1/4-filled band 
with the bond distortion of Fig. 1(e) (short dashed line). 
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j (second chain) 

FIG. 3. The z-z spin correlations between sites 2,3,4,5 on 
the first chain and sites 1 - 12 on the second chain (see text). 
The bars in the top panel show the charge densities Nj on the 
sites of the second chain. 
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FIG. 4. Schematic of the quasi-2d SDW in the correlated 
1/4- filled band. The heights of the arrows have the same 
meaning as in Fig. 1(f). 
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